Radon transform on the cylinder and tomography of a particle on the circle 
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The tomographic probabihty distribution on the phase space (cylinder) related to a circle or 
an interval is introduced. The explicit relations of the tomographic probability densities and the 
probability densities on the phase space for the particle motion on a torus are obtained and the 
relation of the suggested map to the Radon transform on the plane is elucidated. The generalization 
to the case of a multidimensional torus is elaborated and the geometrical meaning of the tomographic 
probability densities as marginal distributions on the helix discussed. 
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I. INTRODUCTION 

The Radon transform [l| is the key mathematical tool 
to reconstruct the tomographic map of both the Wigner 
quasidistribution 0, Q of a quantum state 0, 
and the probability distribution on the phase space of a 
classical particle [1, 0- In the quantum case, this sub- 
ject not only motivated refined theoretical approaches 
based on the maximum likelihood estimation, in order to 
extract the maximum reliable information llOl . but also 
interesting experiments with photonic states [lli. photon 
number distributions l^l and (helium) atoms [13j , focus- 
ing in particular on the reconstruction of the transversal 
motional states. A scheme has been also proposed in or- 
der to obtain the tomographic map associated with the 
longitudinal motion of a neutron wave packet Re- 
cent progress on the quantum aspects has been driven 
by modern experimental techniques and good reviews on 
these topics can be found in [l5j . 

The tomographic map provides the symplectic tomog- 
raphy [TgI of quantum states connected with the sym- 
plectic transform on the phase space (the plane for 
one degree of freedom) and this map can be consid- 
ered as a specific tomographic version of the star-product 
quantization [3 ■ Notice that this interpretation of 
the Radon transform differs from the original motiva- 
tion for the Radon transform in a essential way. The 
genuine Radon transform was introduced as an integral 
transform defined over submanifolds of the configura- 
tion space, more specifically geodesies (i.e., straight lines 
in M^), whereas in symplectic tomography it is rather 
associated to Lagrangian submanifolds of phase space. 
Therefore, although we consider motion, this is instru- 
mental for the identification of the relevant phase space, 
but the actual motions (the solutions of the associated 
Hamilton equations) do not appear in the definition of 



the Radon transform. 

If we consider the classical motion of a particle on a 
circle and its trajectory in phase space (a cylinder of ra- 
dius R) , the motion is described by the time dependence 
of the coordinate q{t) = R4'{t), where 4'{t) is the an- 
gle defining the point on the circle. The angular mo- 
mentum J is the longitudinal coordinate of this motion 
in the phase space. In the presence of fluctuations, the 
particle state is not determined by the two coordinates 
q and p (or (j) and J), but rather by their probability 
distribution function f{q,p) (or f{(l),J)) on the phase 
space. The invertible tomographic map of this distribu- 
tion onto the tomographic probability distribution en- 
ables one to determine the state of the classical particle 
by means of the probability density ujf{X, u), that de- 
pends on a random variable X and two parameters n and 
v. The parameters fx and v label the reference frame in 
the phase space, when the random position X of the par- 
ticle is measured. The reference frame is obtained from 
the initial one by first squeezing the axis q q' — sq, 
p ^ p' ^ s^^p, and then performing the rotation 
q' q" = q' cos 6 + p' sin 9, p' ^ p" ^ ~q' sin 6 -\- p' cos 9 
(see formulae below). Thus the real parameters /i and v 
are expressed in terms of the squeezing s and rotation 9 
as /i = scosd and v = s'mO. The tomographic Radon 
transform maps the probability density, that depends on 
two random variables — position and momentum — onto 
the tomographic probability distribution of only one ran- 
dom variable. 

The case of the motion on the circle can be viewed in 
the limiting case i? — > oo as the motion on the line. Since 
the tomographic map for the classical motion on the line 
is known (and it is very similar to the standard Radon 
transform), it is interesting to address the question of 
whether it is possible to describe the classical motion on 
the circle by an analogous probability density distribu- 
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tion depending on one random variable and some extra 
parameters. The motion that we consider is purely in- 
strumental in order to identify the phase space and does 
not provide us with specific trajectories on which we in- 
tegrate to perform a Radon transform. In fact, not only 
we will discuss the Radon transform of functions depend- 
ing on points on the cylinder (which, to the best of our 
knowledge, has never been presented in the literature), 
but also intend to study how to construct the map of the 
positive probability density distributions living on the 
phase space onto the family of the positive probability 
distributions of random variables living on the helices. 
We will address only the classical motion since the quan- 
tized version of the map, that is known for the motion on 
the line, needs additional consideration for the motion on 
the circle, due to specific properties of compactification 
in one dimension when one goes from the plane to the 
cylinder. The analysis carried out in this paper might be 
therefore very relevant for tomography in quantum me- 
chanics, where we would like to integrate on Lagrangian 
submanifolds to have marginals on the transversal La- 
grangian leaf and therefore it becomes relevant for us to 
understand what is the space of all Lagrangian subman- 
ifolds and the transversal ones. 

The aim of this work is to introduce an invertible tomo- 
graphic map of probability distributions on phase space 
of a particle moving on the circle onto the probability 
marginal distributions on the helix of the cylinder (to- 
mograms). The paper is organized as follows. In Section 
ITTl we review the symplectic tomographic approach for a 
free particle moving on the line. Section IIIII introduces 
the tomographic map for functions on the phase space 
(cylinder) of the free particle moving on the circle. We 
consider an explicit example in Section HVl The multidi- 
mensional generalization is considered in Section |Vl In 
Section IVIl we look at the limit of the tomographic map 
for the particle moving on the torus when the radii of 
the circles tend to infinity and show that in this limit we 
get the symplectic tomographic map corresponding to the 
standard Radon transform. Perspectives and conclusions 
are presented in Section IVIII 



II. SYMPLECTIC TOMOGRAPHY 

Let us consider a function f{q,p) on the phase space 
{q,p) e of a particle moving on the line q E M.. The 
Radon transform as originally formulated solves the fol- 
lowing problem: to reconstruct a function of two vari- 
ables, say f{p,q), if its integrals over arbitrary lines are 
given. 

In the {q,p) plane, a line is given by the equation 

X- nq-vp = 0. (1) 
By using the homogeneity we may write 

X - cos 6'(7 - sin 6*13 = 0. (2) 



Thus, the family of lines has the manifold structure M x §, 
with § the unit circle, X e M and e [0, 27r]. There is an- 
other way to recover this manifold structure which turns 
out to be useful for generalizations to higher dimensions. 
The Euclidean group E(2) acts transitively on the set 
of lines in the plane with a stability group given by the 
translations along the line itself. Therefore the family of 
lines is given by E(2)/M, i.e. R x §. 

The action of R x S may be visualized in the follow- 
ing way: a fiducial line passing through the origin may 
be translated along the normal to the line to generate a 
family of parallel lines. See Fig. [TJ Afterwards, by using 
the rotation group we may rotate this family of parallel 
lines into any other family of parallel lines. As the two 
actions commute, we may also rotate first and then trans- 
late. Thus, we may consider the set of all lines passing 
through the origin and parametrized by the angle and 
then translate each one along the normal. 

It is interesting to observe that 

M2 = E(2)/S, Rx§ = E(2)/R. (3) 

The Radon transform maps TiM?) into JF(M x §), where 
T is a, suitable class of functions that depends on the 
physical setting (for our purposes, is enough). The 
set of lines can be parametrized by two numbers: the 
distance from the origin, d G R, and the angle with re- 
spect to the p = Q axis, 6 e [0, 27r). Any point in R^ can 
be then parametrized by 

(q,p) = (s cos 6,ss\i\9) + (—0? sin 9, d cos 9) , (4) 

where s is the parameter running along the line defined 
by d and See Fig. [T] 

The Radon transform is defined by 

/ + 00 
/(scos6l-dsin6',ssin6l + ficos6')ds. (5) 
-oo 

The inversion formula, as given by Radon, amounts to 
consider first the average value of F on all lines tangent 
to the circle of center P = {q,p) and radius r, namely, 

1 f^^ 

Fp (r) = — / F{qcos9 + p sin 9 + r, 9)d9 (6) 
27r Jo 

and then 

/(g,p) = -- / ^dr. (7) 
TT J r 

The Radon transform maps a (suitable) function on the 
plane into a function on the cylinder. Some conditions 
that guarantee the invertibility and continuity of the map 
were studied by Radon himself [H, John 119|, Hel gason 
[lOl and Strichartz [2lj. 

It is possible to write the Radon transform in the affine 
language (the so-called tomographic map) [l|, [S^ 

u}f{X,^i,v) = {5{X - fiq - lyp)) 

= / fiq,p)S{X - ^q~ iyp)dqdp, (8) 

JR2 
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FIG. 1: Tomography on the plane. 



where S is the Dirac function and the parameters 
X, G M. We notice that 



but also 



fj, V 



1 
1 



-1 

1 



V 



= /.iq + vp. 



(9) 



(10) 



This means that the argument in the Dirac delta function 
may be considered either as a Euclidean product or as a 
symplectic product. Equivalently, one might consider the 
Euclidean or symplectic Fourier transforms. 

Another remark is the following. The full linear inho- 
mogeneous group acts transitively on the family of lines 
on E^. Instead of E(2) as a privileged group, we may 
consider 

SL(2,M) = Sp(2,K) = IGL(2,R)/(R2 x M), (11) 

where SL, Sp and IGL are the special linear, symplectic 
and inhomogeneous linear groups, respectively. The M- 
group in the "denominator" gives dilations while gives 
translations. Because Sp(2, R) is not abelian, it can be 
generated by two types of transformations: rotations 



cos 9 sin 9 
— sin 9 cos 9 



and "squeezing" transformations 



s 




(13) 



The action of the "squeezing" transformation maps lines 
into lines, while preserving the area of the triangle. The 
further action of the rotation group will change the angle 
formed with the q axis. One may show that the Radon 
transform is equivariant with respect to the action of 
SL(2,R) or E(2); both of them preserve the measure on 

R2. 

The inverse transform of ^ reads [l|, 

mp) = [ c.,(^,M,-)e'^^-'^— (14) 

In polar coordinates, fi = rcos9, v = rsin^, the in- 
version formula takes the form of the standard inverse 
Radon transform: 



27r 



dXdd 



Ljf {X, COS e, sin 9)K{9, q, p) , (15) 



f{q,p) - 
with 

K{9,q,p)^sin9 e-"'rcos8-^prsin9^^^ ^^g) 

Jo 

and where we made use of the homogeneity of luj {X, fi, i/) 



ujf{XX,Xp,,Xiy) = —Uf{X,^i,iy), 



(17) 



that is a direct consequence of dH]). If the function f{q,p) 
is a probability density distribution on the phase space 
of a classical particle, i.e. 



fiq.p) > 0, 



f{q,p)dqdp = 1, 



(18) 



(12) 



also the function LUf{X, p,, ly) is nonnegative and is called 
a symplectic tomogram or the "Radon component" of the 
distribution function f{q,p) (analogously to the Fourier 
component of a function). The Radon component con- 
tains the same information on the state of the particle 
evolving on the phase space as the initial distribution 
function. Summarizing: 

w/(X,^,i^) > 0, ujf{X,iJ,,:y)dX = 1, V/x,j/, (19) 
Jm 

the family of tomograms depends on the two real param- 
eters fi and 



III. TOMOGRAPHY ON THE CIRCLE 

In order to extend the preceding tomographic analysis 
to particles confined to compact domains there are two 
alternative definitions, following two different strategies. 
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A. First definition: tomograpiiy on the strip 

Let us choose for definiteness an interval of width 27r. 
The configuration space 



/=[0,27r) 



(20) 



yields the phase space / x R (a strip). To consider this 
case it is convenient to deal with the parametrization of 
lines given by § x E, where M is the translation along the 
normal to the line. If we consider the intersection of the 
lines with the selected strip, it is still possible to consider 
the treatment of the planar situation, where in addition 
the measure dqdp is multiplied by the characteristic func- 
tion of the strip. 

The state of a classical particle moving in the interval 
in the presence of fluctuations is associated with a distri- 
bution function f{q,p) > 0, satisfying the normalization 
condition 



f{q,p)dqdp = 1. 



(21) 



Ixl 



In this case, the symplectic tomogram ([8]) specializes to 



LOf{X,ll,v)= f{q,p)6{X-nq-vp)dqdp, (22) 

J/xR 

with Xjfi,!/ e M. One easily checks nonnegativity and 
normalization like in Eq. 1)19^ : 



LUfiX, n, iy)>0, / ujf{X, v)dX ^ 1, V/x, ly. (23) 
Jr 

The inverse transform, still given by (|14|) . yields a func- 
tion 



.f{q,p) = xi{q)f{q,p) 



(24) 



(xi being the characteristic function) , that vanishes iden- 
tically outside the strip, i.e. f{q,p) = for g ^ /. 

On the other hand, a function / on the strip / x M can 
be extended to a periodic function /27r over the whole 
plane defined by 



= 5^/(<z-27rfc,p)x/+2.fc(9), (25) 



fce2 



where the periodicity, f2Triq + 27r,p) = /2ir('Z,p), is ap- 
parent and we used Eq. ([24|) in the second equality. 

The phase space has become a cylinder § x K, where 
§ = R/(27rZ) is the unit circle. In order to emphasize 
this change of geometry, we will denote the position of 
a particle on the circle by the angle (j) and its angular 
momentum by J. The state of a classical particle moving 
on the circle in the presence of fluctuations is associated 



with the distribution function /((/>, J) = f2Tr{q 
(f>,p — J), satisfying the normalization condition 



d(l)dJf{(f), J) = 1. 



(26) 



Due to the periodicity f{<p + 2kn, J) = f{(j), J) {k £ Z), 
in the inversion formula (jl4p . the Fourier integral over fi 
will be replaced by a Fourier series. Therefore, it follows 
that, in order to reconstruct f{(f>,J), in ((^ only the 
tomograms ujf{X,m,i') with to g Z are really needed. 
Thus, we define 

Lof\x, m, ly) = {d{X - m(j} ~ vj)) 



d<pdJf{,p,J)SiX -m(f)-vJ), (27) 



7x1 



where X, G R and m E Z. In Eq. (P7|) one integrates 
along the family of one-step segments of helices: X = 
m(f> + vJ with < < 27r and X/v — 211111 /v < J < 
Xjv. The choice of this family implies the choice of one 
particular fiber of the cylinder along which each segment 
is discontinuous. In fact, observe that if the 0-domain of 
integration in ([TT)) is changed, say to 



/ + a = [a, 27r -|- a). 



(28) 



one gets different families of tomograms labeled by a 
gauge a. 



Wj"' (X, TO, v) 



dcj)dJf{cj),J), 



/qXI 

which are related to ([27|) by 
_ „^ - ,.,(0) 



(X, TO, v) ~ Lo)^ fi-^ ~ ma, TO, J/), 



(29) 



(30) 



where Taf{q,p) = f{q + Q:,p) is a horizontal translation 
of /. Notice that, due to the periodicity of /, the hori- 
zontal tomogram, with to = 0, is gauge invariant, namely 
tjj"' (X, 0, z^) = uj''j \x,0,iy). Moreover, all families are 
obtained by restricting a G [0, 2tt). In fact one gets 



(31) 



for fc e Z. The gauge a is the anomaly of the chosen fiber 
of the cylinder S x R. See Fig. [2l^a). One easily checks 
nonnegativity and normalization in the form 



W^"^ (X, TO, iy)>0, / W^"' {X, TO, iy)dX = 1, 

Vto, ly, a. 



(32) 



Let us emphasize again that in these formulas, unlike in 
Eqs. (Ull) and m G Z. 

The inverse transform is 



i{X —rrKp—vJ) 



dXdv 



i2ny 



(33) 
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(b) 



a = 2tt + a 



FIG. 2: Tomography on the circle: (a) strip; (b) cylinder. 

Indeed; by making use of the Poisson formula 

e""^ = 27r ^ (5(<?!) - 2m7r) = 27r 62^ {(f)) , (34) 

where 6t is the T-periodic delta function, 



Srit) = S{t (modT)) = 



Sit), 



T = 



(35) 



one gets 



E 



dXdi^ 



e'(^-""^-'''^)4"\x, m, J/) 



(27r)2 
diPdKfiiP^K) J2 



mG2 



xe*(^-™*-''^)(5(X - mV' - ^if) 

^ ^ pim{Tp-(t)) 

d^dKf{i^,K)Y, 



27r 



<5(A" - J) 



dipdKfii^, K)52.{,p - cj>)5{K ^ J) 



la, XR 
/(</', J), 



(36) 



as required. 



B. Second definition: tomography on the cylinder 

When we restrict our attention to periodic functions, 
we are identifying the line at with the line at 27r. In 
this way lines become helices. In this situation, however, 
a new phenomenon takes place: translations along the 
"normal" will map the helix into itself, for translations 
which are integer multiples of 27rtan0 (See Fig. [3]). The 
set of different helices is, therefore, parametrized by an 
angle 9 G (— 7r,0) and the intercept S [0, 27r). Notice 
that the value 6 ~ Q does not correspond to an helix 



J 




^ r((^,e\ 
















231 



FIG. 3: Phase space and relevant variables for the tomogra- 
phy on the strip. 



but to an infinite family of circles "parallel" to the base 
circle. Thus, the set of helices is a trivial bundle with 
fiber § and base manifold S\{0}. where we can use as 
coordinates the slope and intercept (0, 0) or the slope 
and the shift with respect to the helix crossing the origin 
i.e. {e,r[(j),e)) with r{(j),e) = (27r- 0)tan6' . 

Thus, in this setting we would define the Radon trans- 
form as going from functions on [0, 2-k\ x R to functions 
on S X (§\{0}) . It seems clear that only specific ap- 
plications may suggest to use one or the other. For X- 
ray tomography the integration along "segments" may be 
appropriate. For quantum tomography we may want to 
integrate along maximal Lagrangian submanifolds to get 
the marginals along transversal Lagrangian submanifolds 
out of the Wigner function on the full phase space. 

For these reasons wc introduce a different tomographic 
probability distribution: let 

LOf{X, m, v) = ((527r,„(X - mcj) - vj)) 

d4>dJ f{4),J)52TTm{X - m4) - vJ), (37) 

xR 

where X, S R, m e Z and §1 = § is the unit circle. 
Observe that ([57)1 is independent of the (/)-domain of in- 
tegration, due to the periodicity of the integrand. By 
plugging ([55)1 into ([57]) we get for m e Z\{0} (and an 
arbitrary a G R) 



ujf{X,m,iy) 



d(j)dj 



/„xl 



;/(0, J) ^ 5{X -m(t)-vJ - 2TTmk) 
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V / d(j)dj 

fcgZ-^(-fc«+27rfe)xR 

x/(0 - 27rfc, J)5{X -m<j)- vJ) 



d(t)dJf{(t),J)5{X-m(t)-vJ), (38) 

while, for m = 0, 

Cjf{X,Q,u)= ( d(j)dJf{(j),J)S{X-iyJ). (39) 

In conclusion, here we integrate over the whole helix, 
while the previous Eq. ([29|) was integrated on a single step 
of it. Notice also that translations along the line X—m(t)~ 
vJ preserve the measure. By using the homogeneity in 
Eq. p7p we may consider the quantity X/m—cj)— {v /m)J 
that implies 27rm-periodicity of (2; , 

Co f{X + = Co f{X,m,v). (40) 

Therefore, the tomogram lives on a family of cylinders 
labeled by the integer m. See Fig. [2Kb). 
The inverse transform is given by 



(41) 

where §_,„ = §,„ = R/(27r7TiZ) is the circle of radius \m\ 
and §0 = the real line. When / is nonnegative and 
normalized as in (j26p . one easily obtains 

ojf{X,m,i^)>0, / Cdf{X,m,iy)dX = 1, Vm, i/.(42) 

The proof of Eq. (|¥T|) goes as follows 

^ / S^e^^-™^--)^,(X,^,.) 

dHKf{^,K)Y.j ^ 

ixR ^ggJS„,xR l^T'^j 

xe'(^-"'^-''-^)j2™(^ - mV^ - vK) 
d^dKf{ij,K) 
dXdv 



Si xR 

xE 



(27r)2 



Sixl 



^im{tp-(p) 



2tt 



-6{K - J) 



d^dKfi^, K)S2A^ - ^)SiK - J) 

ISi xR 

= /(0,J), (43) 

where we made use of Poisson formula ([34]) and of the 
equality 

f dXe'(''-"'*-''-^'>S2^,n{X-mi'-,yK) 
= V /" dXe*(^-"*-'^^)(5(X - - :/if). (44) 



It is easy to see how the transforms ([^^ and ([57)1 are 
related. Indeed, for m e Z\{0} we get from (|38p 

tD/(X, m, i^) = ct;|"''(X — 2'Kmr, m, :^), (45) 



while, from (139] 



ivf{X,0,,,)^J;'\x,0,,y) 



(46) 



Incidentally, this relation can be used to give an alterna- 
tive proof of the inversion formula (|4ip . In fact, from the 
equality 



/ Cuf{X,m,iy)e'^dX 



J2 I ujf\x ^2TTmr,m,v)e'-''dX 

rgZ la,a+27rm} 



ujf{X,m,i^)e'^dX. 



(47) 



which is trivially valid for m = 0, the inversion formula 
(f33ll translates into diT]) . 



C. A few comments 

A few comments are in order. If the configuration 
space is an interval, the phase space will be a strip and 
a "free" particle bouncing back and forth will move on 
a rectangle. If we impose periodic boundary conditions 
we get circles parallel to the base. Clearly, if we want 
to consider the quantum case, we have to integrate the 
Wigner function on Lagrangian subspaces and get the 
marginals, out of which we should be able to "recon- 
struct" the function. As wc know, we need a "large" fam- 
ily of such marginals, perhaps parametrized by the sym- 
plcctic group, to be able to reconstruct the "state," i.e. 
the original Wigner function 0, i, i, IH Hi] . This view- 
point differs from the original Radon formulation based 
on the set of geodesic lines of the plane M? as Riemannian 
space (for the two dimensional case), whereas in our case 
the relevant lines are the Lagrangian lines of the symplec- 
tic plane as phase space of the one dimensional particle. 
In the Radon case the picture is dynamical while in the 
symplcctic case is purely kinematical. 

In our "classical" setting, we asked a similar question, 
i.e. how to reconstruct a classical distribution function 
on phase space by means of its integrals on a family of 
one-dimensional subspaces. In some sense the fact that 
the family is parametrized by two numbers appears as a 
necessary condition for the reconstruction to be possible. 

Finally, it appears that the two ansatz considered in 
this section yield two different phase spaces. It is rea- 
sonable to expect that what is a suitable function in one 
situation, need not be suitable for the other one. There- 
fore the two proposals may coexist, once it is clear that 
they represent different physical situations. In general. 
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they will yield different results. In a way, physics will 
decide which transform better matches the problem at 
hand. 



IV. GAUSSIAN EXAMPLE 



Let us consider as an illustration the particular case 



1 



(27r)3/2 



(48) 



which is properly normalized, /g^g/ = 1- The Radon 
transform (|27p yields for m 7^ 



— — e ^ 5{X -m<t>-vJ) 



1 



d4) 



(27r)3/2 



exp 



{m(t)- Xf 



47r|m| y/TT 



1 



l=^(a-f +27r) 



%/2| 



r("-lS) 



47r|m| 



VV2|J 



a h 27r 

m 



-erf 



( 1^ 



V%/2|i^| 



(49) 



where erf(x) is the error function. On the other hand, if 
m = 0, 



uf\x,Q,v) 



— _e ^S{X-vJ) 
7„xR (27r)3/2 



1 



(27r)i/2|jy 



X2 

2i^- 



V. TORUS TOMOGRAPHY 

The generalization to many particles is straightfor- 
ward. Let us consider > 1 classical particles, each 
moving on its own circle. The system state is described 
by a probability distribution function /(0, J) > satis- 
fying the normalization condition 



T"xR" 



d(AdJ/(0, J) = 1, 



(52) 



with coordinates = (0i, . . . , 07v) G = (§)^ on the 
A^— torus and angular momenta J — (Ji, . . . , Jn) G M^. 
The tomogram of the torus is defined by 



N 



N 



= f d$dJf{l J) TT 5{Xk - nikA - i^kJk), (53) 

with G and m G Z^. The inverse transform 

reads 



dXdO 



(27r) 



2N 



N 



XLU 



f\x,rn,,7)l[e 



i(Xk-mkipk-VkJk) 



(54) 



k=l 



The tomograms ^'•"^(X, m, t7) and a)(X,TO, j7) are ob- 
tained analogously, as N dimensional generalizations of 
mi and (1371). 



It is easy to verify that the inverse Radon transform ([33]) 
permits to recover the original function (|48p . 

On the other hands, the tomograms along the helices 
read (m 7^ 0) 



Ldf{X,m,v) = 



^dJ 



dJ 



e 2 (52t(X — m0 — i^J) 



(27r)3/2 



(27r)3/- 



■ exp 



(77l(/)- X)^ 
2^^2 



27r|m|' 



(51) 



while, for m = it coincides with ((50|) . (I)/(X, 0, i^) = 
w^"^(X,0, 1/). Note that Eqs. (ggD and dH]) satisfy (|i5)) . 

It is clear from this example that the two transforms 
are different. As we stressed before, both being math- 
ematically legitimate, a choice should be motivated on 
physical grounds. 



VI. LIMIT TO THE STANDARD RADON 
TRANSFORM 

Let us discuss now how the formulas for the Radon 
transform (and its inverse) of a function defined on a 
cylinder tend to those of the standard Radon transform 
of a function defined on the plane in the limit of infinite 
radius of the cylinder. To this end, let us first recall 
how the Fourier series of a periodic function /i? (g) with 

p 7?/2 

period R and normalization /_^y2 ./i?('7)o''Z = 1 becomes 
the Fourier integral when R ^ 00. The Fourier series 
reads 

27r 



and its coefficients are given by 



1 /-^/^ 

^ J~R/2 



^i27Tmq/ R 



dq. 



(55) 



(56) 
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For R ^ oo the Fourier series becomes the Fourier inte- 
gral representation of the function /(q) — Um/^^oo /h(<z) 
defined on the hne. Thus Eq. ([55|) becomes 



R 



f{q) ^ hm ^Afc— Cfc^e 



(57) 



where A/c km+i — km = 2tt/R, and C(fc) = 
hmCfc^i?/27r. On the other hand, Eq. ([56]) takes the 
form 

C{k) = hm — Cfc„ = — / /(g)e'^-«dg. (58) 

Using these well known limiting relations one can get 
the limit of the tomographic map formulae for the par- 
ticle moving on the circle. For definiteness we will look 
at the tomogram (|27p : the procedure is analogous for the 
other tomograms. We first replace Eq. (|?7|) by a formula 
that takes into account the radius R of the circle. Given 
a probability density f{(j), J) > on the cylinder, by in- 
troducing the new variables q = (j)R/2Tr and p = J and 
setting 



2n f2nq 



R-' \ R 

we have the tomogram (P7)) in the form 

^^f\^, Mm, l^) = {^iX - ^Imq - lyp)) 



(59) 



i?,/2 foo 



B.I2 J- 



fii{q,p)S{X - ^irnq-i^p)dqdp, (60) 



where fim = 2T:m/R with a correctly normalized proba- 
bility density 



/ fB{q,p)dqdp=l. 

-R.I2 J -oo 



r-R/2 

The inverse formula (l33l) reads 



,i(X-/i„g-!yp) 



(61) 



(62) 

with A/i = 2-K / R. In the limit R — > oo, we get formulae 
(IH) and (fT4|) and the tomographic map on the circle yields 
the Radon transform on the plane. 



VII. CONCLUSIONS AND PERSPECTIVES 

We have shown that one can map the probability dis- 
tribution density /((/>, J), defined on a cylinder in terms 



of two random variables (position (j) and angular momen- 
tum J) , onto a family of probability distribution densities 
depending on one random variable X, which is a contin- 
uous coordinate on the helix. The family of helices is 
labelled by the integer number m and the real number v. 
The map is obtained by means of the Radon transform 
extended to the case of a cylinder. 

The Radon transform is closely related to the Fourier 
transform. We pointed out an important specific prop- 
erty of the Radon transform, that is valid both for to- 
mographic maps of functions defined on the plane and 
on the cylinder: in contrast to the Fourier transform, for 
which the Fourier component of the probability density 
is not a probability density, the Radon component of the 
probability density (given on the plane or the cylinder) 
is again a probability density and depends on some extra 
parameters. 

We have also straightforwardly extended the Radon 
transform construction to the classical motion on a multi- 
dimensional torus and shown that the tomographic map 
of probability densities on cylinder becomes the tomo- 
graphic map of probability density on the plane. This 
implies that the two corresponding Radon transforms are 
related to each other, in close analogy to the relation be- 
tween Fourier series and Fourier integrals for functions 
on a circle and functions on a line. One difference should 
be stressed though: while in the Fourier case the limit is 
taken in L^, in the Radon case it is (obviously) taken in 
L^. This is apparent in the manipulations of Sec. IVIl 

The quantum extension of the tomographic map for the 
free motion on a circle requires additional investigation, 
due to the well-known ambiguities in the definition of the 
analogues of the conjugate observables angle and angu- 
lar momentum [25j . Similarly, the extension of Radon 
transforms for curved manifolds in the pre sent and re- 
lated contexts deserves additional study [20| . 
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